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Abstract 

Malaria is one of the most common mosquito borne diseases. Temperature is an 
important factor which affects the life cycle of the mosquitoes and transmission dynamics of 
the malaria disease. In the present work, we use SEIR compartmental model for the human 
population and LSEI compartmental model for mosquito population taking temperature 
dependent parameters. Basic reproduction number, R 0 of the model is computed using Next 
Generation Matrix Method. Stability of the disease free equilibrium and the existence of 
the endemic equilibrium point are discussed by basic reproduction number, R Q . Numerical 
results are carried out with different temperature levels. It is observed that temperature 
affects the transmission dynamics of malaria disease significantly. 
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Introduction 

Change of climate mostly affects the transmission dynamics of vector born disease. The 
disease is transmitted by female Anopheles mosquitoes. It is spreading in almost all tropical 
and subtropical region of the world. About 3.2 million people are at risk of malaria. In 2015 
there were 214 million new cases of malaria and 438000 deaths from malaria. There is not 
any effective and safe vaccine against malaria till date. Control e orts of malaria are based 
on the strategies of reducing the mosquitoes, using the anti-malaria drugs and personal 
protection against mosquito bites. But still million of people are not receiving the services 
[!]• 

Climate factors such as temperature, rainfall, humidity, wind and duration of daylight 
strongly affect the ecological and behavioral features of vector borne diseases [2]. In the late 
19 th century average global temperature was increased by about 0.5°C - 0.6°C and by 2100 
it is estimated that global temperature will rise 1.0°C - 3.5°C [3,4]. 

Change of climate mostly impact on the transmission of vector borne diseases. For 
many diseases it is observed that the transmission of disease occurs within the lower end 
range and upper end range of temperature. These range of temperatures are respectively 
considered as 14°C- 18°C and 35°C- 40°C [4]. In the context of malaria, the female adult 
Anopheles mosquitoes feed more frequently in increasing temperature since blood is 
digested more quickly with higher temperature [8]. Githeko et al. (2015) explained that the 
mosquitos' feeds on blood every 4 days at 17°C and with a temperature increased to 25°C it 
feeds on blood every 2 days. The average lifespan of an adult mosquito is about 21 days and 
it is rapidly decreasing in 30°C- 32°C temperature [3]. Temperature levels above 34°C have a 
negative impact on the survival of both vectors and parasites and survival of adults reaches 
zero at around 40°C temperature [5,6]. Githeko et al. (2015) reported that the development 
of a gambiae larva stops when the ambient temperature is below 16°C and the temperature 
below 14°C leads to their deaths. Generally, the availability and productivity of mosquito 
breeding areas increases by rainfall. Although excess rainfall can rush out the Anopheles 
breeding areas. 

Mathematical model is very useful tool to understand the dynamics of malaria disease 
transmission. The mathematical model for malaria disease transmission was studied by 
Ronald Ross in 1911 using SI compartmental model. His work was extended by Macdonald 
[7,8]. Aron and May described the properties of models defined by Ross-Macdonald. 
Anderson and May reviewed the model incorporating various parameters [9]. Chitnis et 
al. analyzed the bifurcation of a malaria model [10]. Mordecai predicted that the optimal 
temperature for malaria transmission occurs at 25°C [11]. In the present paper we consider 
SEIR model for human (host) and LSEI model for mosquito (vector) with temperature 
dependent model parameters. 

Formulation of the Model 

We follow the compartmental malaria model with temperature dependent model 
parameters to study the effect of temperature on malaria disease. The total human 


Article Information 


DOI: 

Article Type: 
Journal Type: 
Volume: 
Manuscript ID: 
Publisher: 


10.31021/acs.20181107 
Research Article 
Open Access 

1 Issue: 2 

ACS-1-107 

Boffin Access Limited 


Received Date: 
Accepted Date: 
Published Date: 


November 02, 2017 
February 02, 2018 
April 04, 2018 


^Corresponding author: 

G Bhuju 

School of Science 

Department of Natural Sciences (Mathematics) 
Kathmandu University, Kavre, Dhulikhel 
Nepal 

E-mail: bhujug@gmail.com 


Citation: Bhuju G, Phaijoo GR, Gurung DB. 
Mathematical Study on Impact of Temperature 
in Malaria Disease Transmission Dynamics. Adv 
Comput Sci. 2018 Apr;l(2):107 

Copyright: © 2018 Bhuju G, et al. This is 
an open-access article distributed under the 
terms of the Creative Commons Attribution 4.0 
international License, which permits unrestricted 
use, distribution, and reproduction in any 
medium, provided the original author and source 
are credited. 


Adv Comput Sci 


Volume: 1.2 







Journal Home: https://www.boffinaccess.com/journals/advances-in-computer-sciences/index.php 


BoPfin Access 


population N h (t) at time t is divided into four epidemiological 
classes: the susceptible class S h (t) (people who are healthy and may 
potentially get infected with malaria parasite), exposed class E h (t) 
(people who are infected with malaria disease but cannot transmit 
to other), infectious class I h (t) (people who are infected with malaria 
disease and able to transmit the disease) and recovered class R(t) 
(people who have recovered from malaria disease). Thus, 

N„(t)= S h (t) + E h (t) + I h (t) + R(t) 

We divide the mosquitoes into two stages: immature and mature 
(i.e. adult), since temperature affects both stages of the mosquitoes. 
The total mosquito population at time t, N m (t) is subdivided into four 
compartments: immature mosquitoes L(t) (egg, larva and pupa stages 
of mosquitoes), susceptible mosquitoes S m (t) (mosquitoes that may 
potentially get infected with malaria parasite), exposed mosquitoes 
E h (t) (mosquitoes infected with malaria that cannot transmit the 
disease) and infectious mosquitoes I h (t) (mosquitoes infected with 
malaria that can transmit the disease). Then at any time t, 

N(t)= L(0 + S (t) + E(t) + 7 (0 

m h h h 

These state variables and their descriptions are presented in 
Table 1. Based on the description of state variables the SEIR-LSEI 
model is given by the following deterministic non-autonomous, 
system of nonlinear differential equation where T = T (t) denotes 
the ambient temperature at time t. 
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Variables 

Description 

S„ 

Number of susceptible human 


Number of exposed humans with infection 

!i 

Number of infectious humans with infection 

R 

Number of humans recovered from infection 

L 

Number of immature vectors 

S 

m 

Number of susceptible vectors 

E 

_m 

Number of exposed vectors 

I 

_m_ 

Number of infectious vectors 


Table 1: State Variables and their descriptions 


wR 



Figure 1: Flow diagram of the model 


The infection rate a h , to susceptible human is temperature dependent, 
A mh is the probability of malaria transmission to a susceptible 
human per bite from an infectious Anopheles mosquito and b(T) 
is the temperature dependent per capita biting rate of mosquitoes. 
Further j h represents the progression rate of exposed individual to 
infectious class. The infected humans acquire temporary immunity 
at the rate of T to join the recovery class. The recovered humans 
lose the immunity at constant rate \j/ and returns to the susceptible 
class. There is natural death rate in each host class and an additional 
disease related death rate K in the infectious class. 

The female Anopheles mosquitoes take rest for a few days after 
obtaining blood meal from human host. These days the mosquitoes 
digest the blood and develop the eggs. This process depends on 
temperature. Typically eggs hatch within 2 3 days but in cold climate 
this may take up 2 to 3 weeks [12]. The temperature dependent 
parameter (T) is the egg deposition rate of mosquitoes from each 
class. We assume that the immature mosquito population is limited 
by carrying capacity k. So that the logistic growth rate for immature 
mosquitoes is represented by 


A b(t)I A b{t)I 

Where CK* = and a,,, = — - - 

N N 

h h 

The total size of vector population and host population are 
related by 

N h {t) = S h (t) + E h (t) + I h (t) + R(t) 

N m ( t ) = L{t) + S m (t) + E m it) + / (0 ( 2 . 2 ) 

In the system (2.1), is the recruitment rate of human. All the 
recruited individuals are assumed to be uninfected new born when 
they join the susceptible class. (Figure 1). When an infectious mosquito 
bite a susceptible human, the human progresses to the exposed class. 


<t>(T) 


1-- 


(S m +E m +I m ) 


The parameters (7 (T) and CD l (T) represent temperature 
dependent maturation rate and mortality rate of immature 
mosquitoes respectively [13,14]. The single compartment L(t) 
represents the three aquatic stage of the mosquitoes egg, larva and 
pupa. When mosquitoes become adult they enter to the susceptible 
compartment S m (t). Susceptible mosquitoes can become infected 
when they bite infectious human and move to the exposed class 
at the rate of temperature dependent parameter a m and A hm is the 
probability that a bite from a susceptible mosquito to an infectious 
human leads to the infection to the mosquito. The temperature 
dependent parameters p m (T ) and CO (T) represent progression rate 
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of exposed mosquitoes and death rate of adult mosquito respectively 

[15] . Since mosquitoes never recover from infection after they 
are infected (because of short life-cycle of mosquitoes), the vector 
population does not include recover class. 

Temperature dependent parameters 

Temperature affects on the rate of development of immature 
and adult mosquitoes and transmission of malaria. Low 
temperature strongly affects on the development of vector 
larvae. The ambient temperature determines the rate of the adult 
mosquitoes feed on human blood. Small increase in temperature 
would probably produce greater mosquitoes densities, higher 
biting rate and more rapid parasite development in the mosquito 

[16] . This relation gives that the biting rate of mosquito increases 
when the temperature increases and temperature levels above 
34°C have negative impact on the survival of mosquitoes [5]. So, 
when the temperature increases above 28°C, the biting rate of 
mosquitoes decreases strictly and it will be zero in around 34°C 
(Figure 2). Mordecai etal. [11] have experimentally suggested the 
biting rate b(T) as a function defined by, 

b = 0.0002037(7 -11.7) V32.4 - 7; (3.1) 

The egg deposition rate for susceptible, exposed and infectious 
mosquitoes is given by [3] 

0(D = -O.153r+8.617-97.7 (3-2) 


The population density of the mosquitoes increase in the range 
of temperature 17°C-28°C and egg deposition rate also increases in 
this range. Then the population density decreases after 28°C (Figure 

3). 

The temperature dependent parameters, progression rate of 
exposed vectors (3 (T), natural mortality rate of immature and adult 
mosquitoes C0j (T) and CO (T) respectively are defined as [17] 
(Table 2). 

p.(T) = -0.00083r + 0.0447 - 0.487 (3.3) 


lo(T) = - log(— 0.000828r 2 + 0.03677 - 0.522 (3.4) 


w.(7) = - 


8.560 + 20.654 


1 + 


19.759 


(3.5) 


f \ 



Figure 3: Probability of egg survival rate of mosquito versus 
temperature. 

\_y 


The progression rate of mosquitoes (3 (T), increases with 
temperature up to 28°C and then it decreases with increasing 
temperature. This result also reveals that there are more infected 
mosquitoes within the temperature range 17°C-28°C (Figure 4). 
Since the number of eggs laid by female mosquito is temperature 
dependent, the rate of life time egg laid B(T) is defined as [17]. 

w) 

co(T) 

That is, the total number of eggs laid by a mosquito is equal to 
the product of number of eggs laid per female per day (EF D) and the 
average adult mosquito life-span, where 

EFD(T ) = -0.153 r + 8.6LT - 97.7 

The probability that an egg survives to become an adult 
mosquito and the development time from egg to adult mosquito are 
represented by PEA(T) and TEA(T) respectively. Both the parameters 
are temperature dependent and defined as [17] 

PEA(T ) = -0.0009247 72 + 0.4537 7 - A.11 

TEA(T ) = --- 

-0.000947 2 +0.0497-0.552 


"\ 




v 


Figure 2: Biting rate of mosquito versus temperature. 

_ J 


Figure 4: Rate of life time egg laid versus temperature. 
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Parameters 

Description 

Dimensions 

p 

Recruitment rate of the humans 

number X day 1 

^mh 

Transmission probability per contact for susceptible humans 

dimensionless 

hm 

Transmission probability per contact for susceptible vectors 

dimensionless 

P h 

Progression rate of exposed humans to infectious humans 

day 1 

P 

Natural mortality rate for humans 

day 1 

n 

Disease induced death rate 

day 1 

¥ 

Rate of loss of immunity of humans 

day 1 

K 

Carrying capacity of immature mosquitoes 

dimensionless 

b(T) 

Temperature dependent biting rate of mosquitoes 

day 1 

<t>ro 

Temperature dependent egg deposition rate 

day 1 

C7(T) 

Temperature dependent maturation rate from eggs to adult mosquitoes 

day 1 

(O (T) 

Temperature dependent natural mortality rate for adult mosquitoes 

day 1 

(T) 

Temperature dependent natural mortality rate for immature mosquitoes 

day 1 

P n m 

Progression rate of exposed mosquitoes 

day 1 

B(T) 

The lifetime number of eggs laid 

dimensionless 

PEA(T) 

The probability that an egg survives to become an adult mosquito 

dimensionless 

TEA(T) 

The development time from egg to adult mosquito 

dimensionless 

EFD(T) 

The number of eggs laid per female per day 

day 1 


Table 2: Parameters and Their Dimensions 


Finally, the temperature dependent maturation rate from eggs to 
adult mosquito is given as [17] 


B{T)PEA{T) 
TEA(T ) 


(3.6) 


Carrying capacity of human 

From the system of equations (2.1) we have 


dNn 

dt 


< p-jiNh 


Jst^oo,N < 

* P 

Hence the carrying capacity of human is — 


(4.1) 


Equilibrium Points and Stability Analysis 

Disease free equilibrium point 

Disease free equilibrium point is a steady state solution of the 
system of equations (2.1) in the absence of infective population. In 
the present case, 


E =1 =R = E =/ 

h h mm 

So, the system of equations (2.1) has a disease free equilibrium 
point f 

= (S\E\I\R\L\E\V)= -S0,0,0,Z,\S\0,0 

* * * m m 


Where 


dr 

dt 


= <KT) 
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S* m (0 — (&(T) + lu(T))L* 


dS* 

— = (j{T)L -uo{T)S * 


Basic Reproduction Number 

The basic reproduction number, R 0 is the expected number of 
secondary infections produced by a single infective on its infectious 
life. According to the values of the basic reproduction number, disease 
can persist with R 0 > 1 and the disease die out when R 0 < 1. 

We compute the basic reproduction number R 0 , associated with 
the disease free equilibrium point; 

{¥ 0, 0, 0, L*, s m *, 0, 0) 

^ j using the next generation method [18]. Using first four 
equations of the above system of equations, a non-negative matrix F 
of the infection terms and the non-singular matrix V of the transition 
terms are given, respectively. 
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Local stability of disease free equilibrium 

Theorem 1: The disease free equilibrium point for the system of 
equation (2.1) is locally asymptotically stable if R 0 < 1 and unstable 
ifR 0 >l. 

Proof: Jacobian matrix of the system of differential equations 
(2.1) at the disease free equilibrium point 


By the Routh-Hurwitz criteria the disease free equilibrium point 
will be asymptotically stable if R 0 < 1, that is det(H.) > 0 for i = 1, 2, 3, 
4. Here, if R 0 < 1, then R 0 2 < 1; And, 
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The characteristic polynomial is: 
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For the stability of the disease free equilibrium, we have to prove 
all the roots of the polynomial (6.1) lie in the left half of the complex 
plane. For the fourth order polynomial (6.1) to have the roots with 
negative real parts, the Routh-Hurwitz criteria are det(H.) > 0 for i = 
1, 2, 3, 4 


= A 


(P + P +P+P\P(P +P +P) + P(P +P) + PP) 

v 1 2 3 l v 2 3 4' 2 v 3 4 3 4' 
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In the above expression, all the negative terms get canceled. So, 
det(HJ is positive. Since all the determinants of Hurwitz matrices are 
positive, the polynomial (6.1) has the roots having negative real parts. 
Hence the disease free equilibrium point is stable if R 0 < 1. On the 
contrary, if R 0 > 1 then P 1 P 2 P 3 P 4 B < 0. So, A 4 < 0 and hence all the roots 
of polynomial (6.1) cannot have negative real parts. It concludes that 
the disease free equilibrium point is unstable if R 0 > 1. 

Theorem 2: (Existence of Endemic Equilibrium Point). The 
endemic equilibrium point of the system exists if R 0 > 1. 

Proof: Solving the system of equations (2.1) the non-zero 
equilibrium point obtained is (S h e ; E h e ; I h e ; R e ; L e ; S m e ; E m e ; I m e ), where 
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The endemic equilibrium point exists if R 0 2 > 1 i.e R 0 > 1, since 

(p 3 p 4 p 5 -^ /; )>°. 


Numerical Results and Discussion 

In the present work, we use SEIR-LSEI epidemic model of the 
malaria disease with temperature dependent parameters. The 
simulations are carried out to explore the effect of temperature in the 
transmission dynamics of the disease. 

The parameter values for the simulation are considered as shown 
in the Table 3. 

All the figures are drawn for T = 14°C, 16°C, 20°C, 28°C and 32°C. 
The following initial values are used for numerical computation, S h (0) 
= 100000, E h (0) = 4000, I h (0) = 18000, R(0) = 18000, L(0) = 30000, 
SJ0) = 130000, EJ0) = 14000, IJ0) = 4000, N h = 140000, k = 40000. 

Figure 5 shows the dynamics of susceptible human population 
in different temperature levels. It shows the positive impact of 
temperature on the transmission of disease. With the increase in 
temperature level from 16°C, more susceptible humans get infected of 
the disease. When infectious mosquito bites the susceptible human, 
it becomes infected and moves from susceptible class to the exposed 
class. Since biting rate of mosquito reaches its optimum value at 
28°C temperature (Figure 2), the susceptible population size become 
minimum in this time. The population starts to increase due to the 
loss of immunity of recovered human and natural birth of the human 
population. 

Figure 6 shows the dynamics of infectious human with different 
temperature levels. The population decreases with time because 
of the malaria induced death, natural death and development of 
immunities in host population. Infectious host population decreases 
in the temperature level below 16°C and above 32°C than in the 28°C. 
When infectious people get immunity they move to recovery class. 

Larva cannot grow in the temperature below 16°C and in the 
temperature below 14°C larva cannot live. Consequently, the number of 
immature mosquitoes sharply decreases to zero in 14°C temperature 
level (Figure 7). It takes long time for immature mosquitoes to develop 


to adult stage and egg deposition rate approaches to zero around 16°C 
temperature (Figure 3), hence, the number of immature mosquitoes 
decreases slowly and approaches to zero after 30 days (Figure 3). 
Temperature is the main factor to affect the development of immature 
mosquitoes. When temperature increases from 16°C, maturation rate 
of immature mosquitoes increases but time of maturation decreases. 
At around 28°C temperature, immature mosquitoes enter to the 
aquatic stages within 2-3 days. So, in this temperature, the number 
of immature mosquito population decreases sharply in initial stage. 
This is because of the increase of egg deposition causing to the 
increase in immature mosquitoes population reaching near to the 
carrying capacity (Figure 7). 

The biting rate of mosquito is very low in the temperature below 
16°C and is optimum in 28°C temperature (Figure 2). So, the number 
of susceptible mosquitoes decreases slowly in 16°C temperature and 
the number decreases rapidly in 28°C because of infection of the 
disease (Figure 8). The population of infectious vectors increases 
initially due to the interaction of susceptible vectors with the 
infectious hosts. Along with the increase in the biting rate, infectious 
mosquito population increases in 20°C and 28°C temperature levels. 

f \ 



Figure 5: Dynamics of susceptible hosts with different 
temperature levels. 

V_ J 


p -\ 



Figure 6: Dynamics of infectious hosts with different 
temperature levels. 
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Parameters 

p 

P h 

¥ 

T 

p 

n 


Kih 

Values 

2.53 

0.071 

0.011 

1/7 

0.00001814 

0.00000638 

0.24 

0.062 


Table 3: Parameter Values 


M \ 



Figure 7: Dynamics of immature vectors with different 
temperature levels. 

v_ J 

C \ 



Figure 8: Dynamics of susceptible vectors with different 
temperature levels. 

v_ ) 

But, more mosquitoes get infected in 28°C temperature than in 20°C 
because of the biting rate depending on the temperature. Afterwards, 
the population starts decreasing due to death of the mosquitoes. The 
infectious mosquito population decreases in the temperature levels 
14°C, 16°C and 32°C as these temperature levels are not favorable for 
the survival of the mosquitoes. 

Conclusion 

Malaria is an infectious disease which is becoming more prevalent 
worldwide and is spreading in new population of the world. Climatic 
factors like temperature play a significant role in affecting vectorial 
capacities of the mosquitoes and then the transmission dynamics of 


the disease. 

(Figure 9) In the present work, we used epidemic compartmental 
vector host model, SEIR-LSEI to observe the impact of temperature in 
the transmission dynamics and prevalence of the disease. Parameters 
of the model are considered temperature dependent. The numerical 
results show that the temperature levels less than 16°C and greater 
than 32°C have negative impact on the survival of the mosquitoes 
and the transmission dynamics of the disease. So, the people living 
in this range of temperature are at low risk of the malaria disease. 
Meanwhile, it is observed from the simulations that large population 
get infected of the disease when the temperature is mild. Thus, mild 
temperature is favorable for the survival of the mosquitoes and 
transmission of the disease as well. The dimensionless number, basic 
reproduction number R 0 which determines whether the disease dies 
out or persists is computed using Next Generation Matrix method. 
The results show that disease free equilibrium of the model is 
locally asymptotically stable when R 0 < 1 and unstable when R 0 > 

1. The number depends on the model parameters which depend on 
temperature. So, the temperature is responsible in determining the 
value of the number. Thus, the temperature has a significant impact 
on the transmission dynamics of the malaria disease. 
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